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On the basis of the representation of the generalized structure of nucleons a new model of the 
hadron interaction at high energies is presented. The new t-dependence of the generalized par- 
ton distributions (GPDs) is obtained from the comparative analysis of different sets of the parton 
distribution functions (PDFs), based on the description of the whole sets of experimental data of 
electromagnetic form factors of the proton and neutron. Taking into account the different mo¬ 
ments of GPDs of the hadron the quantitative descriptions of all existing experimental data of the 
proton-proton and proton-antiproton elastic scattering from ^/s = 9.8 GeV to 8 TeV, including 
the Coulomb range and large momentum transfers up to —t = 15 GeV^, are obtained with a few 
free fitting high energy parameters. The real part of the hadronic elastic scattering amplitude is 
determined only through complex s satisfying the dispersion relations. The negligible contributions 
of the hard Pomeron and the presence of the non-small contributions of the maximal Odderon are 
obtained. The non-dying form of the spin-flip amplitude is examined as well. The structure of the 
Born term and unitarized scattering amplitude are analysed. It is shown that the Black Disk Limit 
for the elastic scatering amplitude is not reached at LHC energies. Predictions for LHC energies are 
made. 
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I. INTRODUCTION 

A topical problem of modern physics of elemen¬ 
tary particles, the exploring of the dynamics of the 
strong interaction processes at high energies, is con¬ 
sidered in the framework of different approaches us¬ 
ing various models of the structure of hadrons and 
the dynamics of their interactions. Different mod¬ 
els for the description of hadron interaction at large 
distances are developed. They are based on the 
general quantum field theory principles (analytic- 
ity, unitarity, and so on). The relativistic models of 
high ener^ scattering based on the quasipotential 
approach occupy an important place among 

them. Here, the hypothesis about the existence of 
the local smooth quasipotential giving an adequate 
description of high energy scattering processes is es¬ 
sential. In the region of small angles of the scat¬ 
tering the eikonal approach can be used as a conse¬ 
quence of the smoothness of the quasipotential Q. 

The smoothness of the quasipotential is related to 
the dynamics of two-particle interactions and means 
that at high energies the hadrons behave as loose 
extended objects with finite dimensions. 

The elastic hadron-hadron scattering plays an im¬ 
portant role in the investigation of the strong in¬ 
teraction. For the description of the interaction at 
small distances there is the exact theory, QCD, but 
for the interaction at large distances, which is the 
basis for the elastic scattering at small angles, the 
calculation in the framework of QCD is impossible 
at present. These two domains are tightly connected 
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with the experimental determination of the parame¬ 
ters of the elastic scattering and are very important 
for the development of the modern strong interac¬ 
tion theory [^. 

Only in the region of small angles the basic prop¬ 
erties of the non-perturbative strong interaction: the 
total cross section, the slope of the diffraction peak 
and the parameter p{s,t) - ratio of the real part to 
the imaginary part of the scattering amplitude, can 
be measured. Their values are connected, on the one 
hand, with the large-scale structure of hadrons and, 
on the other hand, with the first principles which 
lead to the theorems on the behavior of the scatter¬ 
ing amplitudes at asymptotic energies @ . 

There are indeed many different models for the 
descrhjtion of hadron elastic scattering at small an¬ 
gles [ 3 , Q . They lead to the different predictions for 
the structure of the scattering amplitude at asymp¬ 
totic energies, where the diffraction processes can 
display complicated features This concerns es¬ 
pecially the asymptotic unitarity bound connected 
with the Black Disk Limit (BDL) M and the in¬ 
fluence of the saturation regime on the differential 
cross sections m- 

In the Chow-Yang model [T^, E^l it was as¬ 
sumed that the hadron interaction is proportional 
to the overlapping of the matter distribution of the 
hadrons, and Wu and Yang [l^ suggested that the 
matter distribution is proportional to the charge dis¬ 
tribution of the hadron. Many models used the elec¬ 
tromagnetic form factors of the hadron but, in most 
part, they changed its form to describe the experi¬ 
mental data, as was made in the famous Bourrely- 
Soffer-Wu model [ij. The parameters of the ob¬ 
tained form-factor are determined by the fit of the 
differential cross sections. The authors noted that 
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the form factor is ” parameterized like an electromag¬ 
netic form factor, as two poles, and the slowly vary¬ 
ing function reflects the approximate proportional¬ 
ity between the charge density and hadronic matter 
distribution inside a proton.” 

In paper [l^, it was proposed that the hadron 
form factor is proportional to the matter distribu¬ 
tion. The matter distributions in the hadron are 
tightly connected with the energy momentum ten¬ 
sor [^. In [13 , it was noted that ” the gravitational 
form factors are related to the matrix elements of 
the energy-momentum tensor in a hadronic state, 
thus providing the distribution of matter within the 
hadron”. The recent picture of the hadron structure 
is determined by the general parton distributions 
(GPDs) [H, [l^ which include as part the parton 
distribution functions (PDFs). The first and second 
moments of GPDs give two hadron form factors. 

Before the introduction of GPDs a similar repre¬ 
sentation for the scattering amplitude was used in 
work of S. Sanielevici and P. Valin (1984) 

KM^) = \[ (^) 

Here the function L^’^(x) represents the parton dis¬ 
tributions and Tp’^{k) represents the transfer mo¬ 
mentum dependence. Note that the whole structure 
of the amplitude corresponds to the second moment 
of GPDs. 

Usually, models of high energy hadron interac¬ 
tion include a different kind of leading Reggions: 
one or a few pomerons (including the soft and hard 
pomerons), the odderons with intercept equal to the 
pomeron (maximal odderon) or with intercept close 
(or less) unity and sometimes the spin-flip ampli¬ 
tude. The effect of the hard pomeron contribution 
on the elastic differential cross sections is very im¬ 
portant for understanding the properties of QCD in 
the non-perturbative regime [^. Note that p{s,t) 
of the hard pomeron is essentially larger than p{s, t) 
of the soft pomeron. In [l^, it is suggested that 
such a contribution can explain the preliminary re¬ 
sult of the TOTEM Gollaboration (2^ on the elastic 
proton-proton differential cross sections. 

In our high energy general structure (HEGS) 
model [ 2 ^, the real part of the hadronic amplitude 
is determined only through complex s satisfying the 
cross symmetric relation. In the framework of the 
model, the quantitative description of all existing 
experimental data at 52.8 < ^/s < 1960 GeV, in¬ 
cluding the Goulomb range and large momentum 
transfers 0.0008 < |t| < 9.75 GeV^ , is obtained with 
only 3 fitting high energy parameters. The compar¬ 
ison of the predictions of the model at 7 TeV and 
preliminary data of the TOTEM collaboration are 
shown to coincide well. In [ 2 ^, the contribution of 
the hard pomeron in the elastic scattering at small 
angles at high energies was examined. It was found 


that such a contribution is invisible in the existing 
experimental data including the new data of LHC. 
In [ 2 ^, they came to the same result. 

In the framework of the model, only the Born 
term of the scattering amplitude is introduced. Then 
the whole scattering amplitude is obtained as a re¬ 
sult of the unitarization procedure of the hadron 
Born term that is then summed with the Coulomb 
term. The Coulomb-hadron interference phase is 
also taken into account. The essential moment of 
the model is that both parts of the Born term of the 
scattering amplitude have the positive sign, and the 
diffraction structure is determined by the unitariza¬ 
tion procedure. 

Now we present the extended variant of the HEGS 
model , based on the assumption that the hadron 
interaction is sensitive to the generalized parton dis¬ 
tributions (GPDs), whose moments can be repre¬ 
sented in the form of two different distributions: 
charge and matter, separately. Hence, this model 
used the exact electromagnetic and matter form- 
factors determined by one function - generalized par- 
ton distributions (GPDs). Both the form factors are 
independent of the fitting procedure of the differen¬ 
tial cross sections of the elastic hadron scattering. 
Note that the form of GPDs is determined, on the 
one hand, by the deep-inelastic processes and, on the 
other hand, by the measure of the electromagnetic 
form factor from the electron-nucleon elastic scatter¬ 
ing. We support this picture by a good description 
of the experimental data in the Goulomb-hadron in¬ 
terference region and large momentum transfer at 
high energies by one amplitude with a few free pa¬ 
rameters. 

The differential cross sections of nucleon-nucleon 
elastic scattering can be written as the sum of dif¬ 
ferent helicity amplitudes: 

^ +4i$5n.(i) 

The total helicity amplitudes can be written as 
$,(s,t) = F;^(s,t) , where 

comes from the strong interactions, F!f™(s,t) from 
the electromagnetic interactions and <p{s,t) is the in¬ 
terference phase factor between the electromagnetic 
and strong interactions 

The structure of the paper is as follows: first, in 
Section 2 the basis of the first variant of the HEGS 
model is discussed shortly. 

In Section 3, the hadron form factors are analyzed 
by the new form of the general parton distributions 
(GPDs) with taking into account two forms of PDEs, 
which give the best descriptions of the electromag¬ 
netic form factors of the proton and neutron. As 
a result, the new forms of the electromagnetic and 
matter form factors are obtained. 

Section 4 is devoted to the study of all existing 
experimental data of the proton-proton and proton- 
antiproton elastic scattering from yi = 9.8 GeV 
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to 8 TeV, including the Coulomb range and large 
momentum transfers up to —t = 15 GeV^ in the 
framework of the model with a few free fitting high 
energy parameters. 

In section 5, the structure of the elastic hadron 
scattering amplitude, obtained in the framework of 
the model, is analyzed. First, the Born term of the 
scattering amplitude is discussed with its cross-even 
and cross-odd separate parts. Then the obtained 
overlapping function is considered in the impact pa¬ 
rameter representation. And, finally, the full form 
of the scattering amplitude, obtained after integra¬ 
tion over the impact parameter, is considered. Espe¬ 
cially, we examine the energy and momentum trans¬ 
fer dependence of the slope of the scattering ampli¬ 
tude as the Born term and as the full term of the 
scattering amplitude. Finally, the obtained results 
and the comparison with other models and some pre¬ 
dictions of our model are discussed in Section 6. 


II. THE HIGH ENERGY GENERAL 
STRUCTURE (HEGS) MODEL 


F 2 {t) = f dx 

Jo u,d 

In the t-dependence of GPDs in the form 

W{x,t) = q{x)nf exp[a+ t]; (3) 

(1 — x')‘^ 

£‘^{x,t) = qix)sf exp[a_ ——— t]; 

was researched. The function q{x) was chosen at the 
same scale = 1 as in [s^, which is based on the 
MRST2002 global fit [H. 

For ^ = 0 one has the second moment of GPDs 

f dx x'^[H{x,t) zb 8(x,t)] = Ah{t) zb Bh{t). (4) 
^0 


u.d 


The integration of the second moment of GPDs over 
X give the momentum-transfer representation of the 
form factor. It was approximated by the dipole form 

Hence, the Born term of the elastic hadron ampli¬ 
tude can be written as 


The model is based on the representation that 
at high energies the hadron interaction in the non- 
perturbative regime is determined by the reggenized- 
gluon exchange. The cross-even part of this ampli¬ 
tude can have two non-perturbative parts, possible 
standard pomeron - P 29 and the cross-even part of 
three non-perturbative gluons P^9- The interaction 
of these two objects is proportional to two different 
form factors of the hadron. This is the main as¬ 
sumption of the model. Of course, we cannot insist 
on the origin of the second term of the scattering 
amplitude. However, in any case, it has the cross¬ 
even properties, positive sign and the slope is the 
same as for the odderon. The second important as¬ 
sumption is that we chose the slope of the second 
term four times smaller than the slope of the first 
term. All terms have the same intercept. 

The form factors are determined by the general 
parton distributions of the hadron (GPDs) [^. The 
first form factor, corresponding to the first momen¬ 
tum of GPDs is the standard electromagnetic form 
factor - G{t). The second form factor is determined 
by the second momentum of GPDs -A{t). The pa¬ 
rameters and t-dependence of the GPDs are deter¬ 
mined by the standard parton distribution functions, 
so by the experimental data on the deep inelastic 
scattering and by the experimental data for the elec¬ 
tromagnetic form factors (see El). 

The electromagnetic form factors can be repre¬ 
sented as the first moments of GPDs with ^ = 0 



F^^°™(s,t) = hi G\t) Fa{s,t) (I + ri/s°-5) (5) 
+h2 A^{t) Fbis,t) (1-b r2/s‘’-®), 

where Fa{s,t) and Fb{s,t) have the standard Regge 
form 

Fa{s,t) = Fb{s,t) = ‘. (6) 

The slope of the scattering amplitude has the loga¬ 
rithmic dependence on the energy, B{s) = a' ln{s), 
with a' = 0.24 GeV“^ and s = se~™j sq, sq = 
1 GeV^. The final elastic hadron scattering ampli¬ 
tude is obtained after unitarization of the Born term. 
So, first, we have to calculate the eikonal phase 

I (7) 

and then obtain the final hadron scattering ampli¬ 
tude 

Fh[s,t) = is j b Jo{bq) r(s, 6) db (8) 

with 

F(s,6) = I - exp[x(s,6)]. (9) 

The model has only three high-energy fitting pa¬ 
rameters and two low-energy parameters, which re¬ 
flect some small contribution coming from the dif¬ 
ferent low-energy terms. 

We take all existing experimental data in the en¬ 
ergy range 52.8 < Xs < 1960 GeV and the region of 
the momentum transfer 0.0008 < —t < 9.75 GeV^ 
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TABLE I: Experimental data of the electromagnetic 
form factor) 


N points 

111 

196 

87 

Proton 

r^v 

ryp 

References 

[37-43] 

[37,39,40,44-46] 

[38-43] 

[38,39,44,47-49] 


neutron 


13 

Gl 

[50-56] 



[57,58] 

38 

/^n 

[59-63] 

6 


[52,64] 



n 


FIG. 1: The sum of xf of the descriptions of the proton 
and neutron electromagnetic form factors by different 
PDFs over an increasing number of free parameters [Eqs. 
(10) and (11)]. 


of the elastic differential cross sections of proton- 
proton and proton-antiproton data [s^. So we in¬ 
clude the whole Coulomb-hadron interference region 
where the experimental errors are remarkably small. 

As a result, one obtains J2Xi/^ — 1-8, where 
N = 975 is the number of experimental points. Note 
that the parameters of the model are energy inde¬ 
pendent. The energy dependence of the scattering 
amplitude is determined only by the single intercept 
and the logarithmic dependence on s of the slope. 

In the framework of this model the quantita¬ 
tive description of all existing experimental data at 
52.8 < ^/s < 1960 GeV, including the Coulomb 
range and large momentum transfers (0.0008 < \t\ < 
9.75 GeV^), is obtained with only three high-energy 
fitting parameters. Hence, the model is very sensi¬ 
tive to any additional contribution. 



FIG. 2: Proton Dirac form factor multiplied by (the 
hard, dotted, long-dashed, and short-dashed lines cor¬ 
respond to our calculations with PDFs A112, StollerOl, 
Rad04, and Kroll, respectively). 


III. GPDS AND FORM FACTORS OF THE 
NUCLEON 

Further development of the model requires care¬ 
ful analysis of the momentum transfer form of the 
GPDs and a properly chosen form the PDFs. In Ref. 
[s^, an analysis of more than 24 different PDFs was 
performed. We slightly complicated the form of the 
GPDs in comparison with Eq. (jd]) , but it is the sim- 

S 5t one compared to other works (for example. Ref. 

): 


t) 


q{x)lf *; ( 10 ) 

q{x)tf *, 


r{x,t) 


qixrfi ‘, ( 11 ) 

qixffi ) *, 


where qix)'^f = q(a:)“j (1. — xY'^’^'^. 

A complex analysis of the corresponding descrip¬ 
tion of the electromagnetic form factors of the proton 
and neutron using the different PDF sets (24 cases) 
was carried out . These PDFs include the leading- 
order, next-to-leading order, and the next-to-next- 
to-leading order determinations of the parton distri¬ 
bution functions. They used different forms of the x 
dependence of the PDFs. The analysis was carried 
out with different forms of the t dependence of the 
GPDs. The minimum number of free parameters 
was six and maximum was ten. 

To obtain the form factors, we have to integrate 
over X in the whole range 0—1. Hence, the form of 
the X dependence of a PDF affects the form and size 
of the form factor. But the PDF sets are determined 
from the inelastic processes only in some region of 
X, which is only approximated to a; = 0 and x = 
1. Some PDFs have a polynomial form of x with 
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TABLE II: The fitting parameters of the GPDs with flavor dependence 


Model 

tu 

±0.02 

±0.01 

eo 

±0.01 

m 

±0.01 

OLH 

±0.03 

CX.E 

±0.07 

Xo 

±0.002 

Zu 

±0.03 

Zd 

±0.03 

x+4p/yo 

ABKM09 

ABM 12 

0.11 

0.13 

0.2 

0.04 

0.09 

0.13 

0.42 

0.41 

0.45 

0.47 

0.57 

0.60 

0.004 

0.002 

0.67 

0.54 

-1.88 

-2.06 

0.91 

0.89 


a different power. Some others have an exponential 
dependence on x. As a result, the behavior of the 
PDFs, when x —>■ 0 or a; —>■ 1, can influence the form 
of the calculated form factors. 

The sets of experimental data are presented in Ta¬ 
ble I. The sets of data have various corrections and 
different methods which take into account the sys¬ 
tematic errors. So we take into account only the sta¬ 
tistical errors. On the basis of this analysis we calcu¬ 
lated the electromagnetic form factors of the proton 
and neutron (using the isotopic symmetry). Then 
we carried out the fit of these calculations and ob¬ 
tained the parameters of the electromagnetic [G(t)] 
and matter [A(t)] form factors. 

The results of the fitting procedure with differ¬ 
ent numbers of free parameters for the whole set of 
PDFs are presented in Fig. 1. We found that the 
best description was given by the PDFs from Refs. 

. In this case, the increase in the number of 
the free parameters leads to a small decrease in 
This means that the x dependence of the PDFs cor¬ 
responds sufficiently well to the u and d distributions 
in the nucleon to reproduce the electromagnetic form 
factors. Note that these PDFs use a special power 
dependence on x. The most stable results (i.e., a 
minimum dependence on the number of free param¬ 
eters with a minimum of y^) are obtained with the 
PDFs ABKM09 [H and ABM12 [gil (see Table II). 

The obtained form factors for the proton and neu¬ 
tron are shown in Fig. 2 and Fig. 3. The form fac¬ 
tors practically coincide for both PDFs used. In Fig. 
2, our results are compared with the other model 
calculations for Fi{t). It should be noted that the 
experimental data for large t were obtained by the 
Rosenbluth method, and our calculations and the 
calculation in Ref. [s^ differ slightly from the ex¬ 
perimental data at large t, but they practically, co¬ 
incide with each other. The ratio of ^Ge/Gm for 
the proton and neutron cases is presented in Fig. 3. 
Our calculations reproduce the data obtained by the 
polarization method quite well. 

On the basis of our GPDs with the ABM12 PDFs 

0 , 

g„(x) = 4.649903 x“°-288(I. - (12) 

_0.593a;-3.607a;^-|-3.718a;^ 

X , 




FIG. 3: The model description of the ratio of the elec¬ 
tromagnetic form factors for the proton and 

for the neutron 


qd{x) = 3.424394x“°-25®(1. - (13) 

1.122a;-2.984a;^ 

X , 

we calculated the hadron form factors using nu¬ 
merical integration, 

Fi{t) = f dx (14) 

Jo 

3 

3 

and then by fitting these integral results with the 
standard dipole form with some additional parame- 
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FIG. 4: The fit of the form factors of the proton: (a) 
(top), the electromagnetic form factor G{t) [Ea. (ll5ll 1 and 
and (bottom) the matter form factor A{t) [Eg. (1171) 1. The 
circles are the moments of the GPDs (shown only every 
tenth point). 


ters for Fi{t), 

^ 4mp - fit _1_ 

^ Anip — t (1 + q/ai + 9^/02 + g^/a|)2 ’ 

The matter form factor 

A{t) = f X dx (16) 

Jo 

is fitted by the simple dipole form 
= (A2 - t)2 ■ 

The results of the integral calculations and the fit¬ 
ting procedure are shown in Fig. 4. Our description 
is valid up to a large momentum transfer with the 
following parameters: 

ai = 16.7, al = 0.78, a| = 12.5 and A2 = 1.6. 

These form factors will be used in our model of the 
proton-proton and proton-antiproton elastic scatter¬ 
ing. 


IV. EXTENSION OF THE HEGS MODEL 

The obtained form factors differ slightly from 
those used in our previous work [^ . Hence, we 
have to make a new fit of high-energy data, includ¬ 
ing now the new data of the TOTEM Collabora¬ 
tion SI13- As was noted in our previous work, the 
model also describes low-energy data qualitatively. 
Now we include in our fitting procedure additional 
experimental data on the pp and pp elastic scatter¬ 
ing up to 8 TeV > y/s > 9.8 GeV. As a result, the 
amount of experimental data increases by a factor of 
3.5 (from 980 to 3416). This gives us many experi¬ 
mental high-precision data points at small momen¬ 
tum transfer, including the Coulomb-hadron inter¬ 
ference region where the experimental errors are re¬ 
markably small. Hence, we can check our model con¬ 
struction where the real part is determined only by 
the complex representation of s = s/soexp{—i'K/2). 
We do not include the data on the total cross sec¬ 
tions crtot(s) and p(s), as their values were obtained 
from the differential cross sections, especially in the 
Coulomb-hadron interference region. Including such 
data decreases but it would be double counting 
in our opinion. We also do not include the inter¬ 
polated and extrapolated data of Amaldi and 
only include their original experimental data. 

As in the old version of the model, we take into 
account only the statistical errors in the standard fit¬ 
ting procedure. The systematic errors are taken into 
account by the additional normalization coefficient 
which is the same for every row of the experimental 
data. It essentially decreases the space of the pos¬ 
sible form of the scattering amplitude. Of course, 
it is necessary to control the sizes of the normal¬ 
ization coefficients so that they do not introduce an 
additional energy dependence. As we will see later 
(Tables HI and IV), the distribution of the coeffi¬ 
cients has the correct statistical properties and does 
not lead to a visible additional energy dependence. 

Such a simple form of the scattering amplitude in 
the huge region of energy requires careful determina¬ 
tion of the slope of the scattering amplitude. As was 
noted in Ref. 0), analytic S'-matrix theory, per¬ 
turbative quantum chromodynamics, and the data 
require that Regg e trajectories be nonlinear com¬ 
plex functions |69l . [z3- The Pomeron trajectory has 
threshold singularities, the lowest one being due to 
the two-pion exchange, required by the t-channel 
unitarity Ref. 0. This threshold singularity ap¬ 
pears in different forms in various models (see 0 )- 

In the present model, a small additional term is 
introduced into the slope which reflects some possi¬ 
ble small nonlinear properties of the intercept. As a 
result, the slope is taken in the form 

B{s,t) = (ai -\-kqe~^^ *'^)Ln{s). (18) 

This form leads to the standard form of the slope as 
t —>■ 0 and t —>■ oo. Note that our additional term 
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at large energies has a similar form as an additional 
term to the slope coming from tt loop examined in 
Ref. and recently in Ref. [^. The basic Born 
amplitudes were taken in the old form, [Eqs. (5) 
and ([n|)] with fixed ai = 0.24 GeV“^ and A = 0.11. 
Taking into account the Mandelstam region of the 
analyticity scattering amplitude for the 2 —>■ 2 scat¬ 
tering process s + u + t = Arrip we take Sq = 4mp 
where irip is the mass of the proton. 

Then, as we intend to describe sufficiently low en¬ 
ergies, possible Odderon contributions were taken 
into account: 

Fodd{s,t) = ± hodd ^^{t) Fb{s,t), (19) 

where hodd = ihdt/ (1 — r^t). 

Just as we supposed in the previous variant of the 
HEGS model that Fb{s,t) corresponds to the cross¬ 
even part of the three gluon exchange, our Odderon 
contribution is also connected with the matter form 
factor A{t). Our ansatz for the Odderon slightly 
differs from the cross-even part by some kinematic 
function. The form of the Odderon working in all 
t has the same behavior as the cross-even part at 
larger momentum transfer, of course, with different 
signs for proton-proton and proton-antiproton reac¬ 
tions. It has a large preasymptotic part, and as a 
result, such a preasymptotic part of the cross-even 
part is practically not felt. Hence, the Born term of 
the elastic hadron amplitude can now be written as 

= hi Flit) Fais,t) (1 + ri/s°-5) (20) 

+h2 A^it) Fbis,t) 

±hodd A^{t)Fbis,t) (1 -h r2/s°-®), 

where Fais,t) and Fb{s,t) are the same as in the 
previous variant of the model [see Eq. ([6])]. 

The analysis of the hard Pomeron contribution in 
the framework of the model [1^ shows that such 
a contribution is not felt. For the most part, the 
fitting procedure requires a negative additional hard 
Pomeron contribution. We repeat the analysis of 
in the present model and obtain practically the 
same results. Hence, we do not include the hard 
Pomeron in the model. 

At large t our model calculations are extended up 
to —t = 15 GeV^. We added a small contribution 
from the energy-independent part of the spin-flip 
amp litude in a form similar to that proposed in Ref. 

[zl- 

FsfisA) =hsfq^Flit)e-^‘f‘^\ ( 21 ) 

It has two additional free parameters. Of course, 
at lower energy we need to take into account the 
energy-dependent parts of the spin-flip amplitudes. 
However, this requires including additional polariza¬ 
tion data in our examination which essentially com¬ 
plicates the picture. This is beyond the scope of this 
paper. Such a contribution can be made in future 




FIG. 5: da/dt for pp (top) at y/s = 9.8 GeV and pp 
(bottom) at -^5 = 11.3 GeV. 

works. The model is very simple from the viewpoint 
of the number of fitting parameters and functions. 
There are no artificial functions or any cuts which 
bound the separate parts of the amplitude by some 
region of momentum transfer. 

V. ANALYSIS AND RESULTS 

We included 3416 experimental points were in¬ 
cluded in our analysis in the energy region 9.8 GeV 
< -y/s < 8. TeV and in the region of momentum 
transfer 0.000375 < jtj < 15 GeV^. The experimen¬ 
tal data of the proton-proton and proton-antiproton 
elastic scatterin g ar e included in 92 separate rows of 
32 experiments [34| . including recent data from the 
TOTEM Gollaboration at y/s = 8 TeV [i^. The 
whole Goulomb-hadron interference region, where 
the experimental errors are remarkably small, was 
included in our examination of the experimental 
data (see Tables HI and IV). 

In the fitting procedure by FUMILIM program 
[tJ we calculated the minimum in xf related 

to the statistical errors . The systematic errors are 
taken into account by the additional normalization 
coefficient Uk for the k series (the experiment) of 
experimental data, 

2 _ ''^k ElisA) — Tt{s,t) 

^ ~h ’ ' ' 








tively) and (bottom) at \/s = 8 TeV (circles are TOTEM 
data lH] 


where Ti(s,t) are the theory predictions, including 
the hadronic and electromagnetic parts of the scat¬ 
tering amplitude, and nkEi(s,t) are the data points 
allowing a shift by the systematical error of the k 
experiment (see, for example. Refs. lzi,[ 2 i- 

In the region of small momentum transfer the sys¬ 
tematic errors are on the order of 2 — 5%. For the 
most part, the additional normalization is in the re¬ 
gion 0.95 — 1.05. At large momentum transfer the 
order of the systematical errors is 10 — 20%. In this 
case, the additional normalization is situated in the 
region 0.8 — 1.2. Of course, if one sums the system¬ 
atic and statistical errors, decreases but it 

will give some additional space for changing the form 
of the scattering amplitude. 

Our complete fit of 3416 experimental data points 
in the energy range 9.8 < ^/s < 8000 GeV and the 
region of momentum transfer 0.000375 < —t < 
14.75 GeV^ gives = 1-28, with the pa¬ 

rameters hi = 3.67; /i 2 = 1.39; hodd = 0.76; ko = 
0.16; Tq = 3.82, and the low-energy parameters 
hsf = 0.05; ri = 53.7; ra = 4.45. 

Obviously, for such a huge energy region we have a 
very small number of free parameters. We also note 
the good description of the GNI region of momen¬ 
tum transfer in a very wide energy region (approxi¬ 
mately 3 orders of magnitude) with the same slope 
of the scattering amplitude. The differential cross 



FIG. 7: The form and energy dependence of the diffrac¬ 
tion minimum at low energies (long-dashed line: = 

13.4 GeV; thick hard line: ^/s = 18.4 GeV; hard line: 
i/s = 30.4 GeV; dashed line: i/s = 44.7 GeV). 


sections of the proton-proton and proton-antiproton 
elastic scattering at small momentum transfer are 
presented in Fig. 5 at ^/s = 9.8 GeV for pp scat¬ 
tering, and i/s = 11 GeV for pp elastic scattering, 
and in Fig. 6 at ^/s = 7. TeV and ^/s = 8. TeV for 
pp scattering. The model quantitatively reproduces 
the differential cross sections in the whole examined 
energy region in spite of the fact that the size of the 
slope is essentially changing in this region [due to the 
standard Regge behavior log{s] and the real part of 
the scattering amplitude has different behaviors for 
pp and pp. 

The results for the whole energy region for small 
momentum transfer are presented in Table III for 
the proton-proton elastic scattering and in Table 
IV for the proton-antiproton elastic scattering. We 
can see that the values are suitable. We note 
that they are especially small for the high-precision 
FNAL-JINR data which reach a very small size of 
the momentum transfer (up to —t = 0.00037 GeV^ 
[t^ at energies i/s = 13.4, 19.4, 23.4, and 27.4 
GeV. The additional normalization coefficients do 
not show an energy dependence and are distributed 
only statistically. We also include the high-precision 
non-normalized data of the UA4/2 Gollaboration at 
y/s = 541 GeV which reach a very small momentum 
transfer —tmin = 0.000875 GeV^. 

The real part of the scattering amplitude signif¬ 
icantly influences the size and form of the differ¬ 
ential cross sections in the Goulomb-hadron inter¬ 
ference region miii. The second Reggeons also 
have a large slope for the imaginary part and hardly 
change the slope of the differential cross sections. A 
suitable description of both pp and pp experimental 
data in this region supports the determination of the 
real part of the scattering amplitude chosen in the 
model. It should be noted that possible contribu¬ 
tions of the second Reggeons will essentially change 
the form and size of the real part of the scattering 
amplitude. The results presented in Tables III and 
IV show that up to such a low energy we do not feel 
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TABLE III: The proton-proton elastic scattering at small 
t 


y^, 

GeV 

^min 5 

GeV 

^max y 

GeV 

N 


E.xVn 

Uk 

(norm.) 

9.0 

0.00193 

0.04328 

19 

14.4 

0.72 

1.041 

9.3 

0.01268 

0.1147 

28 

21 

0.75 

1.019 

9.8 

0.00115 

0.115 

64 

87.6 

1.35 

1.013 

9.8 

0.0026 

0.12 

23 

31.0 

1.37 

1.074 

9.9 

0.00063 

0.0306 

73 

81.1 

1.11 

1.014 

10.6 

0.00079 

0.01529 

45 

68.0 

1.39 

1.026 

12.3 

0.00066 

0.02928 

58 

46.9 

0.81 

1.018 

13.76 

0.0023 

0.0388 

73 

84.9 

1.16 

1.023 

13.76 

0.035 

0.095 

7 

2.5 

0.36 

1.029 

16.83 

0.0022 

0.0392 

68 

76.9 

1.13 

1.006 

19.42 

0.00066 

0.0315 

69 

79.5 

1.15 

0.996 

19.42 

0.035 

0.095 

7 

12.2 

1.74 

1.008 

19.42 

0.0206 

0.12 

42 

19.9 

0.47 

1.038 

21.7 

0.022 

0.039 

64 

50.1 

0.78 

0.996 

22.2 

0.0005 

0.02978 

64 

55.6 

0.87 

1.007 

23.5 

0.00037 

0.0102 

30 

58.5 

1.95 

1.008 

23.8 

0.0022 

0.0388 

60 

69.1 

1.15 

1.001 

23.9 

0.00066 

0.0316 

66 

76.5 

1.16 

0.988 

27.4 

0.00047 

0.02579 

61 

66.1 

1.08 

0.987 

30.6 

0.016 

0.11 

48 

53.1 

1.10 

1.005 

30.8 

0.0005 

0.0176 

31 

75.7 

2.36 

1.009 

44.7 

0.00099 

0.01856 

40 

51. 

1.16 

1.004 

52.8 

0.00107 

0.05546 

35 

53.2 

1.52 

1.016 

62.3 

0.00543 

0.05122 

23 

31.7 

1.38 

1.005 

7000. 

0.00515 

0.356 

84 

173.4 

2.04 

0.943 

7000. 

0.006 

0.36 

40 

31.4 

0.77 

1.0 

8000. 

0.028 

0.195 

30 

20 

0.7 

0.9 


TABLE IV: The proton-antiproton elastic scattering at 
small t 


^/s, 

^min ? 

^max > 

N 

EyE 

E«xViv 

Uk 

GeV 

GeV 

GeV 




(norm.) 

11.54 

0.0375 

0. 5 

13 

11.5 

0.88 

0.983 

13.76 

0.035 

0.095 

7 

7.4 

1.06 

0.966 

19.42 

0.035 

0.095 

7 

7.3 

1.05 

1.220 

30.4 

0.00067 

0.01561 

28 

28.8 

1.03 

0.974 

52.6 

0.00097 

0.03866 

28 

24.5 

0.875 

0.987 

52.8 

0.0109 

0.0479 

43 

49.9 

1.16 

0.933 

62.3 

0.00632 

0.03821 

43 

55.8 

1.3 

0.996 

541. 

0.000875 

0.11875 

99 

164.7 

1.65 

unnorm. 

546. 

0.00225 

0.03475 

66 

83.7 

1.25 

1.004 

546.6 

0.026 

0.078 

14 

13.86 

1.0 

1.002 

1800. 

0.0339 

0.285 

28 

28.8 

1.03 

1.024 





FIG. 8: da/dt for pp at = 52.8 GeV (top), and pp at 
= 52.8 GeV (midle), and pp at y/s = 7 TeV (bottom) 
[the dashed line (bottom) is the predictions at = 14 
TeV]. 


the essential contributions of the second Reggeons. 
This especially concerns the possible contribution of 
the /o meson. In some models it has an intercept 
essentially above 0.5 and its contribution influences 
the differential cross sections and CTtot(s) and p{s,t) 
in the ISR energy region. Our results practically 
exclude such Reggeons with an intercept above 0.5. 

The form and energy dependence of the diffraction 
minimum are very sensitive to the different parts of 
the scattering amplitude. The change of the sign of 
the imaginary part of the scattering amplitude de¬ 
termines the position of the minimum and its move¬ 
ment with a change in the energy. The real part 
of the scattering amplitude determines the size of 
the dip. Hence, it depends heavily on the Odderon 
contribution. The spin-flip amplitude gives the con- 
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tribution in the differential cross sections additively. 
So the measurement of the form and energy depen¬ 
dence of the diffraction minimum with high preci¬ 
sion is an important task for future experiments. 
In Fig. 7, the description of the diffraction mini¬ 
mum in our model is shown for low energies. The 
HEGS model sufficiently reproduces the energy de¬ 
pendence and form of the diffraction dip. In this 
energy region the diffraction minimum reaches the 
sharpest dip at ^/s = 30 GeV. Note that at this en¬ 
ergy the value of p{s, t = 0) also changes its sign in 
the proton-proton scattering. The pp cross sections 
in the model are obtained by the s —>■ u crossing 
without changing the model parameters. And for 
the proton-antiproton scattering the same situation 
with correlations between the sizes of p(s, t = 0) and 
p(s,tmin) takes place at low energy (approximately 
at PL = 50 GeV). Such a correlation was noted in 
Ref. [12. 

The model reproduces da/dt at very small and 
large t and provides a qualitative description of the 
dip region at —t ~ 1.4 GeV^, for ^/s = 53 GeV 
and for ^/s = 62.1 GeV for the proton-proton and 
proton-antiproton elastic scattering (see top and 
middle panels of Fig.8). The diffraction minimum 
at yi = 7 TeV is reproduced sufficiently well too 
(see Fig.8c). 

In Fig. 9, the description of the differential cross 
sections of the elastic scattering pp at large t and 
different values of s is presented. It is to be noted 
that the calculation of our integrals with complex os¬ 
cillation functions at large momentum transfer is a 
difficult task and requires high-precision of the cal¬ 
culations. In any case, we obtain a quantitatively 
good description of the differential cross sections at 
large t. In this region of t, the contribution of the 
spin-flip amplitude is felt. We take into account 
only the asymptotic part of this amplitude with the 
simplest and energy-independent forms. Although 
it has a small size, its constant is determined suffi¬ 
ciently well, hsf = 0.06 ± 0.004. 

In Fig. 10, the model calculations for the differen¬ 
tial cross sections are shown for the LHC energies. 
Obviously, in the model the difference in the behav¬ 
ior of the differential cross sections between y/s = 7 
and y/s = 14 TeV is not large. For the most part, 
it is reflected in the movement of the position of the 
diffraction minimum to low momentum transfer and 
the increase of the sizes of the differential cross sec¬ 
tions in the minimum and second diffraction maxi¬ 
mum. 


VI. THE STRUCTURE OF THE 
SCATTERING AMPLITUDE 


In the model, only the Born terms of the scatter¬ 
ing amplitude are determined. The separate terms 
and the full Born scattering amplitude have a sim- 



FIG. 9: da/dt at large t for pp at ^/s = 19.4 GeV (hard 
line), y/s = 27.4 GeV, and y/s = 52.8 GeV (dashed line). 
The squares, circles, crossing are the experimental data 
for each case, respectively. 



FIG. 10: The model predictions of da/dt at y/s = 7 
(hard line) and y/s = 14 TeV (dashed line). 


pie form and their imaginary parts do not have any 
oscillating behavior at small momentum transfer. In 
Fig. II and Fig. 12, the parts of the Born terms of 
the scattering amplitude are presented at y/s = 9.8 
GeV and y/s = 7. TeV. At small momentum transfer 
the imaginary part of the Pomeron P 2 g dominates at 
both energies. At small energy, the cross-even and 
cross-odd imaginary parts of P^g are equal to P 2 g 
in the region of momentum transfer 0.7 — 0.9 GeV^, 
and then they dominate. Both parts of Psg have the 
same size at —t > 1. GeV^. But the real part of the 
cross-odd term of Psg dominates at —t > 0.4 GeV^ 
and determines the real part of the full Born term of 
the scattering amplitude. As it has a different sign, 
compared to the real part of P 2 a, the full real part of 
the Born term changes the sign at —t = 0.3 GeV^. 

At y/s = 7 TeV the picture is different. In this 
case, the imaginary part of the cross-odd P^g prac¬ 
tically does not influences the form of the scattering 
amplitude. The imaginary part of the cross-even P^g 
exceeds the imaginary part of P 29 when —t > 0.3 
GeV^ and further determines the full Born term. 

The energy dependence of the imaginary and real 
parts of the full Born term of the scattering ampli- 
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FIG. 11: The magnitude of the Born parts of the pp 
elastic scattering amplitudes at ^/s = 9.8 GeV. (Top) 
The imaginary parts: the sum of all parts (hard line), the 
contribution of (long-dashed line), the contribution 
of the cross-even P3, (dashed line), and the contribution 
of the cross-odd part of P^9 (line with points). (Bottom) 
The same for the real parts of the Born amplitude. 


FIG. 12: The magnitude of the Born parts of the pp 
elastic scattering amplitudes at ^/s = 7 TeV. (Top) The 
imaginary parts: the sum of all parts (hard line), the 
contribution of P2, (long-dashed line), the contribution 
of the cross-even P39 (dashed line), and the contribution 
of the cross-odd part of P^9 (line with points). (Bottom) 
The same for the real parts of the Born amplitude. 


tude is shown in Fig. 13. Note that the imaginary 
part grows with energy at small and decreases at 
large momentum transfer. The real part changes 
the sign and then grows at small momentum trans¬ 
fer but has a small energy dependence at large t. 

The energy dependence of the slope of the full 
Born term of the scattering amplitude is represented 
in Fig. 14. As s is complex, the slope has the real 
and imaginary parts too. The slope changes with s 
most quickly at small momentum transfer. The min¬ 
imum change occurs in the region —t si 0.8 GeV^. In 
Fig. 15, the energy dependence of the ’’basic” slope 
B{s,t) [Eg. dTO)) ] is shown. Its real part has a simple 
form at low energies with a growing maximum at 
high energies that moves towards low t. However, 
the imaginary part has a more complicated energy 
dependence. The maximum of its energy depen¬ 
dence occurs at a momentum transfer larger than 
that for the real part in the region 0.4 < —t < 0.8 
GeV^. The distinction of the form of the slope from 
the constant (linear) slope cannot be explained by 
the absence of the second Reggeons in the model. 
The second Reggeons have real and imaginary parts 
of the same order with a large slope near ai ~ 0.9 
GeV“^. They essentially change the t dependence 


of the differential cross sections at low momentum 
transfer and especially in the GNI region. 

In order to use the unitarization procedure, it is 
necessary to transform the Born term of the scatter¬ 
ing amplitude in the impact-parameter presentation. 
As our Born term has the form factor in a compli¬ 
cated form compared to the simple exponential, this 
transform can be performed only by numerical in¬ 
tegration, [Eq.d?])]. Then, the standard eikonal rep¬ 
resentation is used to obtain the overlapping func¬ 
tion in the impact-parameter representation, [Eqs. 
(0 and [H])] . 

The energy dependence of the corresponding term 
of r(s,6) eq.® is represented in Fig. 16. One can 
see that the black disc limit is not reached at y/s = 7 
TeV [where Im T{s,b = 0) = 0.93] and even at 
y/s = 14 TeV [where Im T{s,b = 0) = 0.95]. The 
real part of r(6, s) (bottom panel of Fig. 16) is small 
and has some important influence only at large im¬ 
pact parameters. Hence, the behavior of the scatter¬ 
ing amplitude will change to the saturation regime 
when the energy grows essentially above the LHG en¬ 
ergies. It is necessary to take this into account when 
the cross sections are extrapolated from the acceler¬ 
ator energies to the energies reached in cosmic-ray 
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FIG. 13: (Top) panel] The imaginary part of the full 
Born amplitude of the the pp elastic scattering ampli¬ 
tudes at ^/s = 7 TeV (solid line), ^/s = 541 GeV (long- 
dashed line), y/s = 52.8 GeV (dashed line), ^/s = 27.4 
GeV (dashed-dotted line), y/s = 9.8 GeV (dotted line). 
(Botton) The same for the real parts of the Born ampli¬ 
tude. 

experiments. 

The corresponding representations for the total, 
elastic and inelastic cross sections are 

(Ttot(s) = 2 J br totXs, b) db, (23) 

with rtot(s, b) = i?e{l - exp[x(s, &)]}, 

tTei(s) = J bTeXs,b) db, (24) 

with rei(s, b) = Re{\ — exp[x(s, 6)]}^, and 

crinel(s) = J brinel{s,b) db, (25) 

with Tineas, b) = Re{l - exp[2x(s, &)]}. 

The energy and impact parameter dependence 
of these values are represented in Fig. 17. The 
rinei(s,6) saturates the unitarity bound up to 6 = 
0.6 fm at y/s = 14 TeV. However, rei(s, 6) and the 
corresponding Ttot = rei-fTinei do not reach the uni¬ 
tarity bound at this energy. In the middle panel of 
Fig. 17, we can see that such saturation takes place 
only at very large energy, y/s = 100 TeV. As a result, 
the maximal growth of AFinei(s 2 , si) = rinei(s 2 ) — 
rinei(s 2 ) occurs at large impact parameters. The 
difference between Arinei(s 2 , Si), Arei(s 2 ,si) and 
AFtot(s 2 ,si) of these values between y/s = 14 TeV 


FIG. 14: (Top) The real part of the slope of the full Born 
amplitude of the the pp elastic scattering amplitudes at 
y/s = 7 TeV (solid line), y/s = 541 GeV (long-dashed 
line), y/s = 52.8 GeV (dashed line), y/s = 27.4 GeV 
(dashed-dotted line), y/s = 9.8 GeV (dotted line). (Bot¬ 
tom) The same for the real parts of the Born amplitude. 




FIG. 15: a) (Top) The real part of B{s,t) fEg. JTSll I of 
the the pp elastic scattering amplitudes at y/s = 7 TeV 
(solid line), y/s = 541 GeV (long-dashed line), y/s = 52.8 
GeV (dashed line), y/s = 27.4 GeV (dashed-dotted line), 
y/s = 9.8 GeV (dotted line), b) (Botton) The same for 
the imaginary parts of B{,t) fEa. (ll8ll 1. 
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FIG. 16: The overlapping function r(s,&) [for the real 
(top) and imaginary (bottom) parts] at ^/s = 9.8 GeV 
(dashed line), = 52.8 GeV (dash-dotted line), = 7 
TeV (long dashed line), y^ = 14 TeV (hard line). 


and y^ = 7 TeV is shown in the bottom panel of Fig. 
17. Though the size of the growth for the elastic and 
inelastic values is very similar, the maximum of the 
growth takes place at different impact parameters. 
The growth of inelastic processes has a peripheral 
character. However, it is due to the saturation of 
such processes in the central region. 

The full elastic scattering amplitude is calculated 
by numerical integration [Eq.(j51)]. The magnitude 
of the real and imaginary parts of the full elastic 
scattering amplitude are represented in Fig. 18 for 
^/s = 9.8 GeV to ^/s = 7 TeV. The imaginary part 
changes by moving its zero from —t = 1.35 GeV^ 
to —t = 0.5 GeV^. The magnitude of the real part 
has zeroes in the examined region of t and mostly 
changes its size, but conserves its form. 

It should be noted that the real part is negative 
at small momentum transfer at y/s = 9.8 GeV. The 
energy dependence of p{s,t = 0) (the ratio of the 
real part to the imaginary part of the scattering am¬ 
plitude) is shown in Fig. 19. Note that we do not 
include the experimental data on p{s,t = 0) and 
(Jtot{s) in the fitting procedure. 

In Fig. 20, the slopes of the hadronic part of the 
full elastic scattering amplitude are represented at 




8 s (GeV) 



FIG. 17: (Top) r(s, b) tot (hard line), r(s, 6)ei(dashed 
line), r(s, b)inei(long dashed line) at y/s = 14 TeV. (Mid¬ 
dle) The energy dependence of r(s, 6 = 0)tot(hard line), 
r(s,6 = 0)ei(dashed line), r(s,6 = 0)inei(long dashed 
line). (Bottom) The differences between r(s, 6)tot(hard 
line), r(s, 6)ei(dashed line), r(s, 6)inei(long dashed line) 
between yfs = 14 TeV - y/s = 7 TeV. 


yfs = 9.8 GeV and yfs = 7 TeV. Obviously, the 
difference between the slopes at small momentum 
transfer is only in the size, though we examined 
such different energies. The nonlinear behavior of 
the slope of the Born term of the scattering ampli¬ 
tude is very weakly reflected in the form of the slope 
of the eikonalized amplitude, but it essentially in¬ 
fluences the fitting procedure. In a larger region of 
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FIG. 18: The magnitude of the pp elastic scattering am¬ 
plitudes after eikonalization of the imaginary part and 
real parts for for ^/s = 9.8 GeV (long-dashed line and 
dash-dotted line, respectively); and for ^/s = 7 TeV 
(hard line and dashed line, respectively). 



Vs (GeV) 


FIG. 19: The energy dependence of p{s, t = 0) (the ratio 
of the real and imaginary parts of the scattering ampli¬ 
tude for pp (solid line) and pp (dashed line) scattering. 
We also show the experimental data for pp (squares) and 
pp (crosses) scattering from Ref. [3^ . 




-t (GeV2) 


the momentum transfer the slope of the differential 
cross section has a significant difference for the low 
and high energies (bottom panel of Fig. 20). This 
is the result of the eikonalization procedure which is 
reflected in the position of the diffraction minimum. 

The comparison of the energy dependence of the 
model calculations of the slope of the differential 
cross sections is represented in Fig. 21. In Ref. [s^, 
the slope was determined as 


f d%b^r{s,b) _ 

^ J debris,b)' 


(26) 


This gives the slope at t = 0, but the experimental 
data for the slope are obtained at small momentum 
transfer and beyond the Coulomb-hadron interfer¬ 
ence region and in some region of t. So, we used the 


FIG. 20: The slope of the full scattering amplitude at 
yV = 9.8 GeV (dashed line) and ^/s = 7 TeV (hard line), 
for the real parts (top), the imaginary parts (middle), 
and the slope of the differential cross sections (bottom). 


standard determination of the slope, 

Bei{s)=Log[^^]/[\h\-M-, (27) 

In the different experimental data, ti and t 2 are 
different. We take some middle points: —ti = 0.04 
GeV^ and —t 2 = 0.05 GeV^. The experimental data 
have large errors. However, the energy dependence 
of our calculations for the most part coincides with 
the energy dependence of the experimental data. 
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Vs (GeV) 

FIG. 21: The energy dependence of the forward elastic 
slope [Eq. (1271) 1 compared to the existing experimental 
data for pp (hard line and open circles) and pp (dashed 
line and squares) scattering. 


TABLE V: The obtained [s^ and predicted sizes of the 
o'tot(s), mb and p{t = 0, s) 


yi, GeV 

^tot —exp 

O'tot 

p{t = 0, s) 

19.42 

38.98 ± 0.4 

39.58 ± 0.8 

-0.005 ± 0.0006 

22.96 

39.42 ± 0.4 

39.89 ± 0.8 

-0.005 ± 0.0006 

52.8 

42.85 ±0.7 

43.15 ±0.5 

0.074 ± 0.005 

541 

62.72 ±0.2 

62.72 ±0.2 

0.128 ±0.005 

1800 

77.3 ± 0.38 

77.3 ± 0.38 

0.127 ±0.02 

7000 

98.0 ±2.6 

97.16 ±0.5 

0.121 

7000 

96.4 ±2. 

97.16 ±0.5 

0.121 

8000 

101. ± 2.1 

99.4 ±0.5 

0.12 

14000 

104 ± 26. 

108.76 ±0.5 

0.1176 

30000 

120. ± 15 

122.7 ± 0.5 

0.11 

57000 

133 ± 23 

135.4 ±0.5 

0.11 


VII. CONCLUSIONS 


We have presented a new model of the hadron- 
hadron interaction at high energies. The model is 
very simple with regards to the number of parame¬ 
ters and functions. There are no artificial functions 
or cuts that bound the separate parts of the ampli¬ 
tude by some region of momentum transfer or en¬ 
ergy. One of the most remarkable properties is that 
the real part of the hadron scattering amplitude is 
determined only by complex energy s that satisfies 
the crossing symmetries. 


The new HEGS model gives a quantitative de¬ 
scription of the elastic nucleon scattering at high 
energy with only five fitting high-energy parameters. 
Our model of the GPDs leads to a good description 
of the proton and neutron electromagnetic form fac¬ 
tors and their elastic scattering simultaneously. A 
successful description of the existing experimental 
data by the model shows that the elastic scatter¬ 
ing is determined by the generalized structure of the 
hadron. The model leads to a coincidence of the 
model calculations with the preliminary data at 8 
TeV. We found that the standard eikonal approxi¬ 
mation [s^ works perfectly well from ^/s = 9 GeV 
up to 8 TeV. The extended variant of the model 
shows the contribution of the ’’maximal” Odderon 
with specific kinematic properties and does not show 
a visible contribution of the hard Pomeron, as in Ref. 


The slope of the differential cross sections at small 
momentum transfer has a small peculiarity and has 
the same properties in the whole examined energy 
region. Such a uniform picture for the slope gives the 
possibility of further research into small peculiarity 
of different cross sections, such as possible oscilla¬ 
tions [s^. Note that we did not see the contribu¬ 
tions of the second Reggeons with a large slope and 


intercept above 0.5 in the examined energy region. 

The obtained value of the total cross sections 
o'tot(s) and parameter p{s, t = 0) are shown in Table 
V. At low energies the model calculation of crtot(s) 
corresponds to the experimental data. The inclusion 
in our fit of the data of the TOTEM Gollaboration 
increases the atot at -v/s = 7 TeV from 95 mb 
to 97 mb, which can be compared with recent data 

iia. 


In Table VI, our model calculations at -v/s = 7 
TeV are compared with experimental data obtained 
at the LHG by different experimental collaborations. 
On the whole, the model calculations correspond 
to the existing experimental data. In our opinion, 
the experimental result of the ATLAS Collaboration 
o'tot(s) is preferable. 

In Table VII, the comparison of our model with 
some others are shown. The first is the old Bourrely- 
Soffer-Wu (-BSWi) model [3. It has a small num¬ 
ber of fitting parameters (seven), but it also takes 
into account the second Reggeon’s contributions. 
We have already noted that in this model the form 
factor is approximated by some function and it rep¬ 
resents the average between the electromagnetic and 
matter form factors. The recent development of this 
model is represented as BSW 2 [s^. In this vari¬ 
ant the number of experimental points increases, but 


TABLE VI: The sizes of the cross sections and the ratio 
of the (Tei/cinei at 7 TeV 



O'tot 


^inel 

-^el/inel 

HEGSo 

95.1 ±0.8 

24.1 ±0.8 

65 ±0.8 

0.37 

HEGSi 

97.2 ±0.5 

24.3 ± 0.5 

66 ±0.5 

0.37 

TOTEM 

98.3 ± 2.9 

24.8 ± 1.3 

73.5 ± 1.9 

0.34 

ATLAS 

95.35 ± 1.36 

24.0 ± 0.6 

71.34 ±0.9 

0.34 
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TABLE VII: The properties of some models ( BSWi BSW 2 [13l, AGN [H, MN [i3|, HEGSo ^,HEGSi 
-present ) 



BSWi 

BSW 2 

AGN 

MN 

HEGSo 

HEGSi 

Nexp 

369 

955 

1728 -t 178 

2654 -f 238 

980 

3416 

^par 

7 + Regge 

11 -1- Regge 

35 

43 

3-t2 

5-k4 

GeV 

24 - 630 

13.4 - 1800 

9.3 - 1800 

5 - 1800 

52 - 1800 

9 - 8000 

At, GeV^ 

0.1 -2.6 

0.1 - 5. 

0.1 -2.6 

0.1 - 15 

8.7 10"'^ - 10 

3.7 10"'‘ 4-15 

{T.X^)/N 

4.45 

1.95 

2.46 

1.23 

1.8 

1.28 


the number of fitting parameters increases as well. 
Both of the models (ours and BSW) used the stan¬ 
dard eikonal unitarization procedure. Another ex¬ 
ample of the model (which does not use the unita¬ 
rization procedure) is based for the most part on Ref. 
[s^ . Its recent variants are presented as AGN[[83] 
and MN [0]. The work is based on many dif¬ 
ferent forms of the parts of the elastic scattering 
amplitude with many additional artificial functions. 
It examined a large number of experimental data 
points but did not include the Coulomb-hadron in¬ 
terference region. The model includes in the fitting 
procedure the experimental data of the total cross 
sections crtot(s) and parameter p[s,t = 0). Such 
an inclusion decreases the total Both models 
summed the statistical and systematic errors and 
have x^/A^ 2.46 and 1.23, respectively. The obtained 
X^ is minimal in the model MN but with a huge 


number of fitting parameters and the inclusion of 
additional artificial functions. Only our model takes 
into account the whole region of momentum transfer 
(3.75 10“^ > |t| > 15 GeV^, which includes the high- 
precision experimental data in the Coulomb-hadron 
interference region. 

A small number of fitting parameters will make it 
possible to explore some fine additional effects like 
possible oscillations of the scattering amplitude and 
find some corrections to the standard eikonal unita¬ 
rization procedure. 
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